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EXPLICIT DETERMINATION OF LATERAL-DIRECTIONAL 


STABILITY AND CONTROL DERIVATIVES BY SIMULTANEOUS 

TIME VECTOR ANALYSIS OF TWO MANEUVERS 

Glenn B. Gilyard 
Flight Research Center 


INTRODUCTION 


The time vector method has been applied to stability and control analysis for more 
than 20 years. As initially developed, the time vector method used a set of known or 
estimated derivatives to investigate aircraft stability (ref. 1). Shortly thereafter, the 
time vector method was applied to the analysis of flight data to determine stability 
derivatives (ref. 2). In past applications of the time vector method, only two unknowns 
in any single equation of motion could be solved for, and therefore only stability deriva- 
tives of the longitudinal short-period mode could be solved for explicitly. Analysis of 
the lateral-directional mode required that the values of some derivatives be assumed 
in order to solve for the remaining derivatives. 

This report extends the time vector method to the simultaneous analysis of two 
maneuvers which differ by some dependent control movement, such as stability augmen- 
tation system feedback. As a result of the simultaneous analysis of the lateral- 
directional mode, stability derivatives and dependent control derivatives can be deter- 
mined explicitly. In the linear case, no estimates of any derivatives are required. 

An example of the application of the technique to flight data is included. 


SYMBOLS 


Physical quantities in this report are given in the International System of Units (SI) 
and parenthetically in U.S. Customary Units. The measurements were taken in Cus- 
tomary Units. Factors relating the two systems are presented in reference 3. 


a^ transverse acceleration at the center of gravity, g 

g acceleration due to gravity, m/sec^ (ft/sec^) 



moments of inertia about X- and Z-body axes, respectively, 
kg-m^ (slug-ft^) 


XZ 


product of inertia referred to the body X- and Z-axes, 

. 2 , i f.2. 

kg-m (slug-ft ) 



Bolling moment t rad/sec ,2 
*X 


L , L , L 
P o P L 1 


constants of L 


VX. 'V 

r o 


constants of L 


0 


N = Yawing moment ? rad/sec 2 


N ,N ,N 
P o P N P Ns 


constants of N 


Ng » N o » Ng 

p o P N 

r o 


constants of N r 


equation 


equation 


equation 


equation 


P,r 

V 


Y = 


Side force 


(Aircraft mass)(V) 


Y , Y > Y 

p o Py P Y. 
r 0 


p o P Y r % 


roll rate and yaw rate, respectively, rad/sec 
velocity, m/sec (ft/sec) 

, l/sec 

constants of Y p equation 


constants of Y„ 


equation 


a 

o 


initial condition angle of attack at the center of gravity, rad 


B angle of sideslip at the center of gravity, rad 

A determinant 


6 


control deflection 


aileron deflection, aileron deflection that produces right roll is 
positive, rad 



phase angle and amplitude ratio of quantity i relative to quan- 
tity j 


<P 


bank angle, rad 



Subscripts: 

P»r,/3,6 


partial derivatives with respect to subscripted variables 


A dot over a quantity denotes the time derivative of that quantity. 


MANEUVER REQUIREMENTS 


Two oscillatory maneuvers which are free of pilot-induced inputs (after initiation 
of the maneuver) and with damping ratios of less than approximately 0. 3 are required. 
The maneuvers must differ by some dependent control variable, for example, roll 
stability augmentation system (SAS) on and roll SAS off. Another possible combination 
would be two maneuvers with the roll SAS on with different feedback gains. 

The dependent control requirement is easily met, since most research and proto- 
type aircraft have independently selectable roll and yaw SAS and some have pilot- 
variable SAS gains . 

The oscillatory mode to be analyzed must be well separated from other mode 
shapes. 


MATHEMATICAL DEVELOPMENT 


The fundamental concepts and the relations that exist for the equations of motion 
in the time vector format are presented in reference 4. 

Since a time vector solution of one maneuver can determine only two unknowns in 
each equation, past applications were explicit only for control-fixed longitudinal short- 
period maneuvers. Analysis of lateral-directional control-fixed maneuvers required 
the assumption of one derivative in each of the equations. However, two properly con- 
ditioned maneuvers can provide the solution of four unknowns, including a control 
derivative, in each equation without any use of estimates. 

The development that follows assumes the general form of a dependent but different 
control input in each maneuver. The same equations apply when one of the maneuvers 
does not have a dependent control input. 

In dimensional form, the linearized lateral-directional equations of motion are 

P --§?*= Y + L pP + L r r + l 6 6 

r - -^p = NJ 3 + N p + N r + N.6 
I z p p r 6 

P = " r + a 0 P + 'f-( < P+ a t ) 

77a = Y J 3 + Y p + Y r + Yc 6 
V t (f p r r 0 
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Analysis of the First Maneuver 


The derivation that follows manipulates the time vector relations of the first 
maneuver to obtain the coefficients of two derivatives in terms of the two remaining 
derivatives. Throughout, the three equations of motion are treated similarly; there- 
fore, only the p equation is used as an example. 

Expressed in time vector format and with p arbitrarily selected as a base, the p 
equation becomes 


P 

P 


4> 


p/P 


*XZ 

*x 


$ r/p L /3 


/ + L 
P/P P 


<t> z + L 

p/p r 


*r/p +L 6 


4 >; 


V P 


( 2 ) 


Since equation (2) represents a two-dimensional vector diagram, the equation can 
be written in cosine and sine components. In matrix notation this is 


^ cos d> n/ 
p 0/p 


sin Vp 


^ COS d> / 

P P/P 


sin Vp 


COS d>. , 7 

P/P I 


XZ 


T r , 

cos <*>• / - L — cos * / - L k 
r/p r p r/p 0 


6 

“ COS A*/ 

P v 6 / p 


sin <j>. 


XZ 


p/p j. 


r T r T 6 

— sin / - li sin d> / — L~ sin / 
p v r/p r p v r/p 6 p 6/p 


(3) 


Solving equations (3) for the derivatives L„ and L as functions of and re- 

r P ^ 

suits in the form 


L„ — L + L n L + L„ Lc 

" r P L. 6 

r 0 

L ~ L + L L +L L r 

P Po P L r P V 6 

r 0 


(4) 


The constant coefficients of equations (4) to be determined from the first maneuver are 


|£j sin ( 4 > 
p|(|P| I 


IfilPI sin (<b , 15 1 sin (^,/p " *r/pj| 


p/p ” *p/p I x |p 


P r 

Sin ( 4 > , - <*> . ) 

P P p/p r/p' 


pIIpI sin ( y P - v P > 


i[-lfl sin< *> 


Wp ’ *P/ P * + ly P S “ Vp Vp’] 


(5) 


leilil 


IpIIpI 8in Vp ' Vp’ 


B 6 

^ - sin (<t , - *. , ) 
p p /3/p 6/p 7 
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where 


A = ^ -I sin (<t> , - ) 

P P| P/P /3/P 

The procedure for determining the derivative coefficients of the r and a^ equations 

is identical to the procedure for the p equation. The coefficients for the yawing moment 
and lateral-acceleration derivatives are then 
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By using the time vector identities discussed in reference 4, only the time his- 
tories of the variables p, r, a^., and 6 are required. 


Analysis of the Second Maneuver 

Because of the dependent control input of the second maneuver, the aircraft response 
is different than in the first maneuver. However, substituting the results of the analysis 
of the first maneuver into the equations of motion leaves only two unknowns in each 
equation. The time vector analysis of the second maneuver thus provides for the ex- 
plicit determination of four unknown derivatives in each equation. 

Substituting the derivative expressions of equations (4) into the p equation (eq. (1)) 
results in the following: 


p -X (v\v y* 


.\j}+/L + 

V p ° 

b v 


L L+L L K \p+Lr + L R 5 
P L r P L 6) r 6 

r 6 b 


(7) 


~L~ 

P 


Expanding and rearranging, 

P - = L /X + L p o P + L r ^ L /3 l ^ + L p L P + + L 6 (X/ + L P Lg P + ^ (8) 

The p equation is now a function of only two unknowns, L r and L^; therefore, a 
vector analysis yields an explicit solution for these derivatives. 

Writing equation (8) in cosine and sine components and matrix notation results in 


L |&! 

X l pi 

r 


cos 4 > o/ + L | P |cos <f , +|-|cos 4> , \ l^lcos <t> , + L l^lcos <I> / + — cos 4>. , 

(J/p P L I p| p/p |p| r/pl /? L |p| Pip P L |p| p/p |p I 6/p 
r / \ 6 6 / 


L_ 1^1 sin 4> , + L l-lsin* , +1 — 1 sin 1^1 sin * , +L 1^1 sin 4> / + 1^1 sin , 

Pj I p | P/p P l IpI P/P I P l r/ Pj i 0 L I Pi P/p P L |P| P/P |P| 6/p 

r r / \ 6 6 


cos <J> 


XZ 


p/p U 


I 


sin <t> 


XZ 


p/p Iv 


cos , - L„ 

r/p P Q P 


cos , - L — cos <J> , 
P/p P 0 P P/P 


_ B T P . 

sin <J>- / - L. ^ sin <J> _ , - L - sin <|> > 

r/p j8 o p /3/p p Q p P/P 


(9) 
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Matrix (9) can then be solved for and L^. Substituting and L~ into 

equations (4), the values of L„ and L are determined. 

r P 

All the derivatives of the r and a^. equations can be determined in a similar manner. 

APPLICATION OF THE TECHNIQUE 


The simultaneous time vector derivative identification technique described in this 
report has been successfully applied to time histories of maneuvers performed with 
three aircraft. 

Figures 1(a) and 1(b) compare flight data from rudder pulse maneuvers with time 
histories calculated from the explicitly determined time vector derivatives . The 
maneuver shown in figure 1(a) was performed with the roll SAS off, and the maneuver 
shown in figure 1(b) was performed with the roll SAS on. 

In the analysis of the roll-SAS-off maneuver, the coefficients of the p and p 
derivatives were determined in terms of the r and 6 derivatives, as shown in 
equations (4). 

Analyzing the roll-SAS-on maneuver with the results of the roll-SAS-off maneuver 
explicitly defines the complete set of stability and control derivatives: 


II 

-5.45 1/sec 2 

Ng = 4.40 l/sec 2 

Yg = -0. 20 

1/ sec 

L = 

-0.78 l/sec 

N = 0. 02 1/ sec 

o 

ii 


P 

P 

P 


L = 

0.92 l/sec 

N =-0.12 1/ sec 

Y =0.01 


r 

r 

r 


L 6 = 

-6.50 l/ sec 2 

N 5 =0.24 l/sec 2 

Y* = -0. 02 

1/ sec 

a 


a 

a 



DESIRABLE FEATURES OF THE TECHNIQUE 


Some of the desirable features of the simultaneous time vector technique are: 

Simultaneous analysis of two maneuvers determines explicitly three stability 
derivatives and one dependent control derivative in each equation. 

It requires only manual analysis of maneuver time histories. 

The derivative expressions can be easily programed on a digital computer, elimin- 
ating the use of actual vector diagrams. 
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The technique can be expanded to use more than two maneuvers, making it possible 
to determine a correspondingly larger number of control derivatives. 

No estimates of any derivatives are required. 

Zero shift biases of the variables are not required. 

A minimum number of variables are required, such as p, r, a^., and 6 for the 
lateral-directional mode. 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., November 10, 1972 
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(a) Roll SAS off. 


Figure 1. Comparison of lateral-directional flight data with time histories 
calculated from derivatives determined with the time vector method. 
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